First passage times of pulling-assisted DNA unzipping 
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We investigate the voltage-driven transport of hybridized DNA through membrane channels. As 
membrane channels are typically too narrow to accommodate hybridized DNA, the dehybridization 
of the DNA is the critical rate limiting step in the transport process. Using a two-dimensional 
stochastic model, we show that the dehybridization process proceeds by two distinct mechanisms; 
thermal denaturation in the limit of low driving voltage, and direct stripping in the high to moderate 
voltage regime. Additionally, we investigate the effects of introducing non-homologous defects into 
the DNA strand. 

PACS numbers: 



I. INTRODUCTION 

Over the last two decades, the rapid development of 
single molecule manipulation techniques has produced 
dramatic advances in chemistry and molecular physics. 
The fields of DNA biotechnology and nanofabrication in 
particular, have benefited from these advances. In turn, 
successes in the field of DNA biotechnology have further 
motivated the study of DNA manipulation techniques, 
with a strong focus on stud ying mechanisms of induc - 
ing DNA dehybr idization f |Lubenskv fc Nelson 2002l| . 
ICocco et al. 2003J ). 

In this study, we consider dehybridizing DNA by 
pulling DNA molecules through a transmembrane chan- 
nel. This work is motivate d in part by t he experimen- 
tal results of N akane et al. |Nakane 2004J . and Bates el 
al. [Bates etal. 2003J. In the experiments of Nakane, a 
probe DNA that extends 14 base pairs beyond the open- 
ing of a trans-membrane channel is used to capture a 14 
base pair target strand (see Fig^. In order to prevent 
complete passage of the probe DNA through the mem- 
brane channel, the probe was attached to an Avidin an- 
chor protein by a 50-base poly-A tail. The target strands 
used in the experiments were either completely homolo- 
gous to the probe strand, or had one nonhomologous "de- 
fect" base. After the probe strand had captured a tar- 
get strand, a transmembrane voltage was applied to pull 
the probe strand through an a-hemolysin channel em- 
bedded in a lipid membrane. As the a-hemolysin channel 
could only accommodate single- stranded Dl>i A molecules, 
any target strand hybridized with the probe DNA would 
need to detach completely before the probe could be fully 
drawn into the channel. As the transmembrane voltage 
was applied, the ionic current through the a-hemolysin 
channel was simultaneously monitored. When any part of 
the probe strand was inside the channel, the channel was 
blocked, and the ionic current fell below the measurement 
threshold. When the probe strand was completely pulled 
through the channel (after the target strand detached), 
ionic currents resumed and were measured. Thus, in the 
experiment of Nakane, the distribution of first passage 



times, T, for the escape of the probe DNA from the a- 
hemolysin channel was measured. 

In order to understand the first passage times observed 
in these experiments, we have produced estimates for the 
thermally averaged mean first passage time (roq) using 
simple two-dimensional and one-dimensional stochastic 
models. We model the DNA hybridization energetics 
using both a simple free energy model where each hy- 
bridized base contributes an equal amount to the to- 
tal free energy, and using free energie s produced by the 
MFOLD 2-state hybridization server |Zuker 2003| . Our 
predicted mean first passage times show the same quali- 
tative features observed by Nakane et al. Under the influ- 
ence of small transmembrane voltages, the first passage 
times are found to depend strongly on the presence and 
energy of defects in the target DNA strand. As the trans- 
membrane voltage is increased, a distinct roll-over in the 
predicted first passage times is observed, and the first 
passage times become relatively insensitive to the pres- 
ence of defects on the target DNA strand. We propose 
that this roll-over in Toq represents a transition between a 
thermally-dominated dehybridization mechanism in the 
low voltage regime, and a driven stripping mechanism in 
the high voltage regime. 



II. MODEL 

Motivated by the do uble-stranded DNA unzipping 
models of Coc co et. al. [Cocco et al. 20031 and Poland 
and Scheraga [Poland fc Scheraea 1970l |. we model the 
process of DNA probe extraction using a two-sided zip- 
per model described in Fig. ^ This model generates 
- — '^^^ — ^^ distinct partial dehybridization states labeled 
(711,77.2), were ni is the number of bases at the 5' end of 
the probe that are dehybridized and have been pulled 
into the transmembrane region, and 712 is the number 
of dehybridized bases at the 3' end of the DNA probe. 
Using this nomenclature, the DNA probe is fully dehy- 
bridized when 711+712 = N, and fully extracted only when 
ni = N and 712 = 0. Note that his model neglects the 
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(ni + 1, 712), and /zi is a phenomenological attempt 
frequency. We note that we could express /xi as 
the product of an attempt frequency and a free en- 
ergy barrier /ii -^ ^e^^'^ 1^'^ . However for the 5' 
bond-breaking transition, as well as all other tran- 
sitions in our model, we pull any dependence on the 
barrier energy into the effective attempt frequency. 
As this state transition involves a net translation 
of a partially charged nucleotide across the mem- 
brane, a drop in the interaction energy between the 
apphed electric field and the probe DNA occurs: 



AF = V(n^ -F 1, 77,2) - T^(ni, na) == qV. 



(3) 



where q is the mean partial charge of a DNA base. 

•5' Bond Formation In the reverse of the 5' bond 
breaking transition, the probe DNA translates one 
base in the 3' direction, and corresponding bases on 
the probe and the target DNA strands hybridize. 
This shifts the system from the (ni, 712) state to the 
(n\ — 1, 712) state, and occurs at rate 



FIG. 1: An example of the DNA detachment experiment. 
The detachment of the length A'^ target DNA is assumed to 
occur via two coordinates n\ and 712, representing unzipping 
from the two ends of the probe DNA. Here, an A^ = 14 target 
strand is shown with a single defect. The probe DNA is shown 
in blue, hybridized target DNA in green, the defect base in 
red, dehybridized bases in brown, and extracted bases in light 
blue. The probability of bubbles forming is ignored. 
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•3' Bond Breaking The hydrogen bonds between two 
hybridized bases located at the 3' end of the probe 
break, thus increasing 71,2 by one. This transition is 
assumed to occur with rate 



effects of internal dchybridization bubbles, which will be 
discussed later in this section. Defining Pini,n2)i't) to be 
the probability of being in state (tt-i, 772) at time t, given 
that the system was initially in state {n[, n'j), we find the 
system evolves according to a master equation 



dP{t) 
dt 



= MP(i). (1) 

In order to produce the elements Mj^ of the transition 
matrix, we define the following five transition mecha- 
nisms between states (71,1,772): 

•5' Bond Breaking The hydrogen bonds between 
bases at the 5' end of the probe DNA break, 
the probe DNA base moves into the transmem- 
brane region, and the state shifts from (711,712) to 
(txi -|- 1,772). This transition path includes a bond 
breaking event, immediately followed by transla- 
tion of the probe DNA, and occurs at rate 



fci = ^le 



-I3{AG+AV) 



(2) 



Where AG is the change in free energy between 
the current state (711,772) and the target state 



Ai2e 



'/3(AG) 



(5) 



•3' Bond Formation Corresponding dehybridized 

bases at the 3' end of the probe strand and the 
5' end of the target DNA strand form a bond, 
decreasing 772 by one. This transition is assumed 
to occur at rate 



r2 = V2 



(6) 



• Free Translation In the case where the DNA probe 

is completely dehybridized (tii + 772 = N), but 
not fully extracted (771 7^ N), the DNA probe 
will translate under the influence of the applied 
potential. The translation changes the state from 
(771, 712) to (771 -|- 1, 772 — 1) or (771 — 1, 772 + 1) and 
occurs with rate 



fc3 = ^3e-^(^^) 
fc4 = A*3e-^(^^) 



(7) 



The two rates k^ and fc4 describe translation in the 
5' and 3' directions respectively. The only signif- 
icant difference between k^ and fc4 is the sign of 



the change in the interaction energy between the 
probe DNA and the appUed electric field. In the 
case of 5' translation AV < 0, while in the case of 
3' translation AF > 0. 

Having defined the possible transitions in the model, 
we arc left with five unknown attempt frequencies 
(/ii, 1^1, /i2, 1^2, Ma); 3'iid an unknown value for the mean 
partial charge per DNA base q. Based on the results 
of Nakane, we set q = 0.4e jNakane 2004J, were e is 
the electron charge. Since free translation and 5' bond 
breaking both involve sliding of the probe DNA strand, 
the time scale for these transitions is set approximately 
by the time scale for translation of the single-stranded 
probe DNA. To find this time scale, we built a simpli- 
fied model of single-stranded DNA translation consisting 
of only the free translation moves described above. Us- 
ing this model, we simulated the complete passage of a 
60-base single-stranded DNA segment through a chan- 
nel 30 DNA bases in length. We then compared the 
mean first passage time for complete transit of single- 
stranded DNA computed from thi s model, with the ex- 
perimental results of Bates et al [Bates et al. 2003J for 
DNA transport through a-hemolysin channels at low 
voltages (1010/xs at 20mV, SSO/xs at 40mV). The best 
match between the computed and experimental results 
occurred when the attempt frequency was approximately 
/^i = 9 X 10^s~^. To determine a value for V2, we fol- 
lowed Cocco et al. |Cocco et al. 2003^ and set V2 to 
equal the inverse self diffusion time for a single nucleotide 
(« 5x lO^s""'^). With rough order of magnitude estimates 
for /xi and 1^2 we then tuned the parameters of our model 
to give a reasonable m atch to the experimental data of 
Nakane [Nakane 2004J for the variation of escape time 
with voltage. After hand tuning, the parameters were 
^1 = 9 X 105s-\ /i2 = 1.4 X 10'^s-\ vi = 3.5 x lO^s-i, 
1/2 = 6 X l{fs-^, and /X3 = 9 x lO^s"^ 

While easy to analyze, our model fails to capture 
a number of aspects of the denaturation and extrac- 
tion process. For example, by assuming the short 
probe-target DNA complex denatures only from the 
ends, we have ignored the effects of bubble formation. 
By bubbles we refer to dehybridized regions of DNA 
that are bounded by at least one hybridized base on 
both sides. While clearly an approximation, we do 
not believe that bubble formation will qualitatively af- 
fect the process of driven DNA hybridization described 
here. Under physiologic conditions, structural fluctu- 
ations in double stranded DNA lead to the formation 
of bubbles that a re typically on the ord er of several 
tens of base pairs [Hanke fc Metzler 2003J . In this pa- 
per we consider the dehybridization and extraction of 
DNA strands only 10-20 bases in length, and thus con- 
siderably smaller than the typical bubble size found 
in DNA at physi ologic conditions. Addi tionally, re- 
cent experiments [Altan-Bonnet el al. 2003f have stud- 
ied the formation of bubbles in small, atypically bubble- 
prone, poly-AT DNA chains approximately 18 bases in 
length. Here, the formation of DNA bubbles approxi- 



mately 2 — 10 bases in length were observed in the DNA 
strand. However, throu gh fluorescence-correlation spec- 
troscopy Altan-Bonnet [Altan-Bonnet el al. 2003f were 
able to measure the typical lifetime of these small bubbles 
to be approximately 50//s. This lifetime is typically much 
smaller than the mean time for the dehybridization and 
extraction of small DNA strands through a-hemolysin at 
moderate driving voltages. 

Additionally we do not model partial-registry binding, 
where the target DNA binds to the probe DNA with an 
offset so that the target DNA "over-hangs" the probe. 
If we assume that the target-probe DNA complex is well 
equilibrated prior to the start of the extraction processes, 
then such out of registry binding is unlikely due to its 
high energetic cost. Moreover, should such binding occur, 
its effects on the DNA dehybridization process would not 
be qualitatively different from the effects of starting in 
a (71.1,^2) state with a similar number of unhybridized 
bases. 

Finally, we do not model the full extraction of the 
probe DNA [i.e., full removal from the a-hemolysin chan- 
nel). As re-binding of the target DNA is unlikely once 
the probe is fully drawn into the channel ( typical tar- 
get concentrations are lO/iM [Nakane 2004J ). complete 
removal of the probe DNA depends primarily on the de- 
tails of the DNA- channel interactions. Previous studies 
[Bates et al. 2001^ , have shown that translation of single- 
stranded DNA 60 bases in length under moderate driving 
(0.G4V) in an a-hemolysin channel occurs on a time scale 
of 50G//S. Rather than model the interactions between the 
probe DNA and channel, we can simply add 500/xs to the 
first passage times for probe escape we compute using our 
model. As we shall see, at low to moderate membrane 
voltages this additional time is small compared to the 
time required to denature the probe DNA. 



III. RESULTS AND DISCUSSION 

We use Eqn^to compute the thermally averaged time 
for the DNA probe to first reach the fully extracted state 
(ni = iV, 712 = 0). To do this, we make the (iV, 0) state 
completely absorbing, and slightly modify Eqn^by set- 
ting P{N,o) (t) = for all t: 



dP, _ V 



Pj{t)M,^j + P,{t)M,^, , i^{N + l). 



(8) 

Note that the first term in EqnJHldoes not include con- 
tributions from transitions out of state {N, 0) while the 
second term does include transitions into state (A^, 0). 
Upon defining S'(t|(ni, nj)) to be the probability that 
the system reaches the absorbing state between time t 
and time t -\- dt, given that it started in state (n'l, n'2) we 
find 
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With S{t\{n[,n'2)) computed from the solution to Eqn 
IHl we can determine the mean first passage time to the 
extracted state 
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FinaUy, to find the thermaUy averaged mean first passage 
time, Toq, we sum the results of Eon II Ul over all possible 
initial states of the probe DNA, 
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where G{ni,n2) is the free energy of state (ni,n2) rel- 
ative to the free energy of the completely dehybridized 
system. The validity of Eqn 1111 is contingent on the as- 
sumption of local thermal equilibrium between the probe 
and the fluid bath. However, as the probe is typically 
maintained in the inserted state for several tenths of a 
second prior to removal JNakane 2004], we expect this 
assumption to hold to reasonable accuracy. 

To investigate how nonhomologous defects influence 
the mean first passage time, we first analyze the process 
of DNA extraction using a simple model where each non- 
defective hybridized base pair lowered the free energy of 
the probe-target complex by a constant amount. A sin- 
gle defective base pair with a binding free-energy lower 
than all other bases was then introduced into the target 
strand. As an example, we considered a short strand of 
eleven bases, and used the double zipper model to gener- 
ate the mean escape time curves displayed in the inset of 
Fig[21 There are two features worth noting; the mean es- 
cape time clearly depends on the location of a defective 
base as well as the energy of the defect, and the posi- 
tion and energy of the defect influence the escape time 
less as the applied voltage is increased. At first glance 
the dependence of the escape time on defect location ap- 
pears counter-intuitive given that the defect energy did 
not vary with the location of the defect along the target 
strand. Nonetheless, the behavior displayed in Fig. |21 
where a defect placed in the center of the DNA strand 
produces maximal reduction in mean escape time, is ex- 
tremely general and rather insensitive to changes in the 
parameter values of the model. 

To understand the source of this effect, consider a DNA 
chain with a defect at position i, when the DNA chain is 
in the {i — l,j) or {j,i — 1) state with j < {N — i -I- 1). 
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FIG. 2: Mean first passage times generated using free en- 
ergies provided by the MFOLD program (see Table U . At 
low voltages, energetic effects account for the majority of the 
difference in escape times between the various defect-carrying 
strands. INSET: Thermally averaged mean escape times of 
a DNA segment 5 bases long. Each hybridized base lowers 
the free energy by 2kT, except for a single defect site (located 
at position i) which lowers the free energy by fcT/10 when 
hybridized. Placing the defect in the center of the segment 
lowers the mean escape time significantly more than placing 
a defect at either end of the segment. 



If the defect is at site i, the probability of the adjacent 
DNA base rehybridizing goes as r/{r + k). Here k is 
the rate at which the defect base dehybridizes, and r is 
the rate of rehybridization of the base adjacent to the 
defect. The effect of introducing a defect is to increase 
k while keeping r constant, reducing the likelihood of 
hybridization of the adjacent base. In the case where the 
defect results in particularly weak binding, the increase 
in k can be quite large, making it highly unlikely that 
DNA will rehybridize beyond the defect base i once the 
defect itself has dehybridized. This shortens the effective 
length of the DNA by min(i — 1,N — i + I) sites. 

We note however that introducing bubbles into the 
model would minimize the variation in escape time with 
defect position. In the presence of bubbles, bases well 
separated from the defect could rehybridize without the 
defect base itself rehybridizing. This would dramatically 
reduce the variation in the mean escape time associated 
with the defect position. Nonetheless, the mean escape 
time of GC-rich strands which are relatively resistant to 
bubble formation may still be sensitive to the position of 
defects in the target strand. 

Turning to a more realistic free energy model, we used 
the the MFOLD 2-state hybridization server (Zuker 20Q3 
to generate free energy surfaces fo r four distinct probe se- 
quences, each 14 bases in length [Nakane 2004J (see Ta- 
ble m . The form of the lifetime verses voltage behavior 
produced by Eqn 1111 for these sequences is displayed in 
Fig. [3 The four voltage response curves clearly separate 



Name 


Sequence 


Energy 


14pc 


5'-GGTGGTTGGTTTGGTT-3' 


-37.7kT 


lA 


5'-GGTGGTTGGTTTGGTA-3' 


-S6.9kT 


IOC 


5'-GGTGCTTGGTTTGGTT-3' 


-27AkT 


7C 


5'-GGTGGTTCGTTTGGTT-3' 


-27.1kT 



TABLE I: DNA target sequences |Nakane 2001 along with 
the hybridization energies predicted by MFOLD (T = 293). 
The probe DNA strand is perfectly complementary to the 
14pc sequence. The difference in the free energies between the 
various defect strands is due to the presence of different defect 
nucleotides (A vs C), and differe nces in the b ases neighboring 
the defects. For more detail see IZuker 2C 



dchybridizcd {i.e., in state (i, A'^ — i)) at time i, assuming 
that no rehybridization events occur between t and t*. 
Rates kl and fcj are defined by Eqns |21 and [51 respec- 
tively. 

Since the probabilities produced by EqnHcontain con- 
tributions from paths through the (ni,n2) state space 
that include rehybridization events, we cannot use the 
solution to Eqn^to compute Qi. However we can com- 
pute Qi by considering the master equation governing the 
transitions between fully dehybridized states (ni, N—ni). 
Defining W(i_Ar_i)(i; {j,N — j)) to be the probabiHty of 
being in the {i,N — i) state at time t, given that the 
system was in state {j, N ~ j) at time i = 0, we have 



into two groups; the 14pc sequence (perfectly comple- 
mentary to the DNA probe) and the lA sequence (with 
a single non-complementary base at the 5' end of the se- 
quence) belong to a long lived set, while the 7C and IOC 
sequences (both with a single non-complimentary base in 
the center of the sequence) form a group with a dramat- 
ically shorter lifetime. We note that the approximately 
four orders of magnitude difference between the defect- 
free and lA escape times and the 7C and IOC escape 
times at zero voltage, is consistent with the free energy 
differences between the various probes. From the stand- 
point of potential applications to DNA sequencing, the 
reduction in free energy produced by a defective base is 
clearly the dominant influence on the mean escape time. 

Returning to the DNA sequences used to produce Fig 
121 we investigated the consequences on the mean probe 
lifetime of varying the 3' detachment rate, and found that 
this rate has a significantly greater influence on the mean 
escape time in the low voltage regime than in the high 
voltage regime. This leads us to propose that the field- 
induced dehybridization can occur through one of two 
mechanisms. In the low voltage regime, thermal denatu- 
ration of the target DNA from the probe is the dominant 
mode of separation. In this region of the lifetime- voltage 
curve, thermally-induced dehybridization from the 3' end 
is at least as significant as voltage driven-bond stripping 
at the 5' end. As the applied voltage is increased, the sys- 
tem enters a transition region where voltage-driven bond 
stripping becomes the statistically favored mechanism of 
detachment. Finally, in the high voltage region stripping 
becomes the primary means of separation. 

To confirm this hypothesis, we find the probability that 
base i (measured from the 5' end of the probe DNA), is 
the last base to dehybridize prior to extraction of the 
DNA probe. We first define g{i.{N-i)-i){t) = 9i{t) to 
be the probability density that the last dehybridization 
event occurs at base i between time t and t + dt. Making 
use of the solution to EquQ] giit) can be expressed as 

g,{t) = PS)k'i^Q^{t,t*) + P,{t)kfQ,+^{tX)- (12) 

Here Qi{t,t*) is the probability the probe DNA is fully 
drawn into the channel by time t* given that it was fully 
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(13) 
In Eqn ^1 there are two distinct sets of absorbing 
states; the first set contains only the fully extracted state 
(iV, 0), the second set contains the states with only one 
hybridized base {i^N — i — 1). As we are seeking the 
probability of reaching the fully extracted state with- 
out having any base rchybridize. we are only concerned 
with the [N, 0) absorbing state. The probability dis- 
tribution for escape from the set of fully dehybridized 
states {i,N — i) to the fully extracted state (A^, 0) is 
Sq{t) = fcr"'V(w„i.i)(i). Then Q,(i,r) is 



Q^{tX) = /o* 'S,{t)dt 



= fcf"'Vo* 'W(N~i,i){t;{hN-i))dt. 

(14) 
Substituting Eqn ^^ into Eqn ^1 and letting t* — > oo, 
we have 
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P,(t)fcf ^fcf -^^ /„°° W^(Ar_i,i)(t; (*, N - l))dt 

+P,(i)r«fcf -^) /„°° W^M-i,i){t; {^. N - i))dt. 

(15) 
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FIG. 3: The last bond probabilities for the 10c probe at var- 
ious voltages. At low voltages the distribution is determined 
primarily by the free energy surface indicating that thermally- 
driven dehybridization is the dominant dehybridization mech- 
anism. With increasing voltage, the effects of the free energy 
surface are rapidly obscured and stripping is the dominant 
denaturation mechanism. INSET: The last bond probabili- 
ties for an 11 base DNA with a defect at position i — 5. The 
configurational free energies were produced assuming a simple 
free energy model with a uniform free energy per bond, and 
uniform attempt frequencies of 10^s~^. 



Note that in the Umit of low target DNA concen- 
trations, once the probe and target DNA have com- 
pletely dehybridized. rchybridization is unlikely during 
the time required to extract the probe DNA. Thus, 
in the case of low target concentrations, the integral 
^3^"^^ /o°° W^N-i,i){t; ih N-i))dt « 1 and EquUilyields 
giit) = Pii,N~i-i){t){r\ + k\). With gi{t), we can com- 
pute the probability, 7^, that the last base to dehybridize 
is base i from the expression 7^ = L gi{t)dt. The ther- 
mally averaged probability that base i is the last to de- 
hybridize is computed from 
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see a similar effect where the free energy surface domi- 
nates the last bond-probability at low voltages. At higher 
voltages the effects of the free energy surface are obscured 
and a stripping mechanism dominates. This is consistent 
with Fig|2 where we see that the mean first passage times 
are most strongly influenced by sequence defects in the 
low voltage limit. 

This transition between a thermally dominated denat- 
uration mechanism and a direct stripping mechanism can 
explain the decreased influence of defects on the mean 
extraction time at high membrane voltages. To demon- 
strate this, we consider the one-dimensional reduction 
of the two-dimensional model of DNA dehybridization 
already presented. Retaining only the 5' hybridization 
transitions, the Master equation for the one-dimensional 
model is 



dP(ni,0) 
Jt 



-{ki{ni) + ri{ni)[l-5n,fl))P{ni,Q)+ 

ri(ni + 1)(I - 5m,N-i)P[ni + 1,0), 
< ni < iV-1. 



(17) 
The one-dimensional model of DNA denaturation de- 
scribed by Egn llTl is equivalent to a one dimensional ran- 
dom walk with a reflecting boundary at ni = and an 
absorbing boundary at ni = N . The first passage time 
for extraction in to the transmembrane channel is then 
jPurv fc Caceresj 
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Referring to the inset in Fig|3| we see that using the 
simple free energy model where each hybridized base low- 
ers the free energy of the system by 2kT produces a 7, 
that is symmetric about the mid-point of the DNA lat- 
tice. This is exactly what is expected in the case where 
all the attempt frequencies in the model are equal and 
thermal denaturation dominates the dehybridization pro- 
cess. As the voltage increases, the behavior of 7, is also 
predictable, with the bases at the 3' end of probe DNA 
becoming more likely to be the last to separate. This 
is consistent with a voltage-driven denaturation mecha- 
nism where stripping of the target DNA occurs. Using 
the MFOLD free-energies for the 10c DNA segment, we 



and the thermally averaged escape time is 
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To make contact with the two dimensional model, we take 
the high voltage limit of Eqn ^| Note that in the high 
voltage limit we expect the extraction times predicted 
by the one and two dimensional models to agree to high 
accuracy. Additionally, we focus on r(0) as this is by far 
the dominant term in Ean ll9l In the limit of high voltage 
ri(ni)/fci(ni) = e'3((G("i+i)-G(«i))-2gy) ^ and 
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rehybridization past the defect base pair, thus shortening 
the effective length of the probe-target complex. 



Here AG(ni) = G{ni + 1,0) - 0(711,0). EqnOOl shows 
that in the high voltage limit, the mean extraction time 
is closely approximated by the time required to indepen- 
dently dehybridize N DNA bases without significant re- 
hybridization occuring. In this limit, a defect only influ- 
ences the extraction process once, when the defect base 
pair is first dehybridized. Conversely, if we consider the 
zero voltage limit. Eon 1181 vields 
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where the second term includes contributions from re- 
hybridization events. The effect of a defect is to reduce 
the ratio ri{j)/ki{j) = {vi/ ^i)e^'^^'^'^^'^\^ and thus at low 
voltages a defect accelerates the extraction process each 
time the defect base pair is denatured. At low voltages 
the base pair is typically denatured multiple times prior 
to probe exaction, and hence defects have a more pro- 
nounced effect on the mean extraction time at low volt- 
ages than they do at higher voltages. We also note that 
the magnitude of Eqn |^ depends strongly on the po- 
sition of the defect in the DNA sequence. As was the 
case in the two-dimensional model, the mean extraction 
time predicted by the one-dimensional model is a min- 
imum when the defect is located in the middle of the 
target strand. As before, we ascribe this position de- 
pendence to a defect's ability to reduce the likelihood of 



IV. CONCLUSION 



We have shown that the dehybridization process crit- 
ical for the voltage-driven transport of DNA through 
membrane channels can proceed through two distinct 
mechanisms. In the low voltage regime, thermal dehy- 
bridization is the dominant mechanism and the details 
of the hybridized DNA's free energy surface set the time 
scale for transport. In the moderate to high voltage 
regime, direct stripping of the complementary DNA off 
the probe DNA is the dominant dehybridization mecha- 
nism. Driven primarily by the applied voltage, this de- 
hybridization mechanism is comparatively insensitive to 
the shape of the hybridized DNA's free energy surface. 
This behavior may have implications for the use of DNA 
sequencing devices based on nanopores. Specifically, an 
attempt to increase the sequencing rate by increasing the 
voltage driving the double stranded DNA through the 
nanopore, may obscure the distinction between different 
DNA sequences. This may limit the maximum rate at 
which sequencing can be performed. 
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